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УДК 681.324 
В.В. МКРТИЧЯН 


КОМПЬЮТЕРНЫЕ МОДЕЛИ СПИСОЧНЫХ ДЕКОДЕРОВ ГУРУ- 
СВАМИ-СУДАНА ДЛЯ ОБОБЩЕННЫХ КОДОВ 
РИДА-СОЛОМОНА И КОНКАТЕНИРОВАННЫХ КОДОВ 


Решена задача разработки компьютерных моделей списочных декодеров Гурусва- 
ми-Судана для обобщенных кодов Рида-Соломона и конкатенированных обобщен- 
ных кодов Рида-Соломона с кодами Адамара: разработан точный алгоритм списоч- 
ного декодирования конкатенированных кодов, получены структурные схемы и 
программные реализации декодеров. 

Ключевые слова: обобщенные коды Рида-Соломона, конкатенированные обоб- 
щенные коды, списочное декодирование 


1. Введение и постановка задачи. Крупным прорывом в теории поме- 
хоустойчивого кодирования было создание М. Суданом в 1997 году принци- 
пиального списочного декодера для кодов Рида-Соломона (РС-кодов) [1]. 
Декодер использует интерполяцию и факторизацию многочленов двух 
переменных над расширением базового поля Галуа и способен с полиноми- 
альной сложностью работать за пределами минимального кодового рассто- 
яния. В работе [2] на основе декодера Судана был получен декодер Гуру- 
свами-Судана для обобщенных кодов Рида-Соломона (ОРС-кодов), имею- 
щий лучшие корректирующие способности. В работе [3, с.177] был предло- 
жен неформализованный метод списочного декодирования Гурусвами-Су- 
дана для ОРС-кодов, специальным образом конкатенированных с кодами 
Адамара (КОРСА-кодов). Как доказано в [4] списочное декодирование Гуру- 
свами-Судана для ОРС-кодов и КОРСА-кодов можно применить для защиты 
тиражируемой цифровой продукции от несанкционированного распростра- 
нения. 

В работе [5] программно реализован декодер Судана для РС-кодов. 
Цель данной статьи — разработка компьютерных моделей списочных деко- 
деров как для более общих ОРС-кодов, так и для КОРСА-кодов. Особенно- 
стью представленной разработки является применение в модели списочно- 
го декодера для ОРС-кодов длиной И и размерностью А алгоритма фак- 
торизации Рота-Руккенштейн [6], позволяющего эффективно, со сложно- 


стью О((/хК + 1ю=а)"1ю=° (г/АК)), проводить факторизацию полиномов 


двух переменных над полем Галуа /”,. 
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2. Необходимые сведения об алгоритме списочного декодирова- 
ния Гурусвами-Судана для обобщенных кодов Рида-Соломона. 


Пусть А”, -— поле Галуа мощностью Ч; Ё’,[х] - кольцо полиномов 


переменной Х над полем /”,; Е’, [х, у] - кольцо полиномов двух пере- 


менных Х и У над полем /”,; В“ '[х] с В [х] - пространство по- 


линомов степени не выше А -1; /А”,” - пространство векторов размерно- 


стью И над полем /",; 4(м,у) - метрика Хемминга, и,уе А.’ [71]; за- 
пись /|Р далее означает, что Л’ делит Р’нацело, где Л. ре Ё.[х, У]. 
Пусть @&@,.--..%, - фиксированное упорядочение элементов /”, ; 
У, ,....У, - фиксированные элементы у (Е №) < а. ОРС-код дли- 


ной И, размерностью А ((/, А) -ОРС-код) можно определить как множе- 
ство векторов (у, р(“,)..... у. р(@&,)), где ВР пробегает 


-—1 

множество информационных полиномов А”, [х]. 
Алгоритм списочного декодирования Гурусвами-Судана ОРС-кодов 
[2] включает два основных шага: шаг интерполяции, на котором по полу- 
ченному слову строится полином двух переменных специального вида, и 
шаг факторизации, где данный полином разлагается на сомножители, по 
которым можно построить список. Входными параметрами декодера яв- 
ляются параметры ОРС-кода: мощность поля Ч, длина Г и размерность 
А кода и некоторый управляющий параметр #(<= | Л-@-= 1) + 1..3) 
При  декодировании на вход алгоритма подается слово 


У= (№,...У,)Е Е,’ в виде сетки {(@, у, );...;(@,,у,)}. Декодер 
производит поиск всех кодовых слов в сфере с центром № радиусом г: 1. 
Выходом алгоритма является список всех информационных полиномов 
ЛО Ех], удовлетворяющих условию: |# | Л(х,) = у > 2. Из [2] 
вытекает, что этот список содержит истинное информационное сообщение. 


Приведем алгоритм списочного декодирования Гурусвами-Судана 
для ОРС-кодов в удобном для нас виде 


АЛГОРИТМ 1./* Вход: 4.7”.А,#; сетка {(@, у, ),....(@,,У,)}. 
Выход: список Л(х).*/ 
Шаг 0. Вычислить параметры: 

т = |+ Ее? + 4(Ё - Ак) (2 - ®)) |+ Ти /= #2 - 1. 





Шаг 1. (Интерполяция) Найти любой полином @(х, у)е ЁЕ.,[х, у], 


в виде 
[2/к | 1- 


С(х,у) 5 у ее ' (1) 


12=0 
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для которого выполняются следующие условия: 
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1. 


ВЮК я 


У в САС АЖ "У * = 0.0) 

2. @(х,у)# 0. 
Шаг 2. (Факторизация) Разложить С’(х, у) на неприводимые сомножите- 
ли. 
Шаг 3. Выдать список всех полиномов Л(Сх)е Ё,[х], таких, что 
(»- /(>х)) является делителем С(х,у), причем Л(х,)= у,, по 
крайней мере в Г значениях #= {1;...; 7}. 

В [2] также имеется "весовая" версия алгоритма 1, которая, полу- 
чая на вход параметры (/’,А)-ОРС-кода, управляющий параметр 


(Е { ку. ий +47 и вектор весов м’ = (м ,..., м’, ) для коорди- 


нат входной сетки {(%,, у, )....,(@,„.У,)}, находит все информационные 


полиномы Л(х), удовлетворяющие условию у влор-» М” 2 Е. Входная 


сетка может иметь длину больше И ‚ включая элементы вида: (ау), 
(а у) ‚ где уя у ‚ что вместе с весами позволяет учитывать вероятно- 
сти появления букв У вточке @ и строить мягкие декодеры и декодеры 
для конкатенированных кодов. Алгоритм модифицируется следующим об- 
разом: на первом шаге параметр Т заменяется на величину 


т, = т] т (т, где Ее {1;..37}, м = тах. М’. Эту вер- 
сию далее будем называть алгоритмом 1”. 

3. Разработка алгоритма списочного декодирования для КОРСА-ко- 
дов. В [3, с. 177], изложен метод списочного декодирования Гурусвами-Су- 
дана для КОРСА-кодов, однако точного алгоритма декодирования не при- 
водится. В этом разделе построен формализованный алгоритм декодирова- 
НИЯ. 

3.1. Специальное конкатенирование ОРС-кодов с кодами Адамара. 


Кодирование КОРСА-кодов. Пусть Р - простое, 7 - натуральное, 
2,>--->2 „ -— фиксированное упорядочение элементов ^”,". Код Адамара 
над полем /”, с инициализирующим параметром 77 задается кодирую- 


ЩИМ отображением 


Ш 
р" Руа = (кал жк ат >), 
т р р 
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т 


и имеет минимальное кодовое расстояние р” - р”"' [7]. Далее этот код 


будем обозначать как (р” , т) -А-код. 
Для описания специального конкатенирования ОРС-кодов с кодами 
Адамара введем ряд обозначений. Пусть Р’ - простое, 7 - натуральное, 


ГЕ {р”;2р";3р";...;р”"}, КЕ {т2т,3т,...тт/р"\, (3) 
И„ - биективное отображение, сопоставляющее элементу Ё,” элемент 
поля /”„ в соответствии с полиномиальным представлением поля, 

Ю=К/т, БЕГ/р”. (4) 
Рассмотрим биективное отображение: 
а: БУ 


= (0) (1) (0-1) Ко-1 
Ди (а) = и „(а )+и,„(а Ее (а . )х° 1 
где а = басы) Г а” = о ‚Е {0;...; Ао ый п} . 
Очевидно, что отображение Х в. определяется формулой: 


Е Е: 
дя (РО)= (и, (Фо (ФО, (Фь)), 


№1 
где р(х)= ро+ Р1х+ ...+ Рь-1х ° _. Рассмотрим отображение: 


ии В; (а) =У „(и (а)), 
где \/ „ - кодирующее отображение (р”,т) -А-кода. Пусть ©, ›---> @ „„ 


— фиксированное упорядочение элементов ^”.». 


Р 


КОРСА-код над полем /”’, , получаемый специальным конкатениро- 


ванием (7.А. ) -ОРС-кода над полем /”„ и (р”,т) -А-кода над полем 


Р 
`, , имеет инициализирующие параметры 77,А,и (см. (3), (4)) и задает- 
ся кодирующим отображением: 


та вы (р. > | „р, 


Ко-1 К 
где Р.(х)Е Ё» [х] - представление сообщения ае РЁ, в виде поли- 


нома над полем /”„, служащее для кодирования "внешним" (7,,Ао )- 


и 
ОРС-кодом над полем ^”,„: р.(х)= Х„„(а); И» - кодирующее отоб- 
ражение "внутреннего" (р” ‚т) -А-кода. 


КОРСА-код размерностью А длиной И, обозначаемый далее как 
(г, А) -КОРСА-код имеет минимальное расстояние 


132 


Раздел “Управление, вычислительная техника и информатика ” 








Я = (1-1/р")(1- (® - П/к)" [3]. Отметим, что в определении КОРСА- 


кода содержится и метод кодирования. 
3.2. Декодирование КОРСА-кодов. Построенный ниже списочный деко- 
дер для КОРСА-кодов состоит из двух основных элементов: внешнего и 
внутреннего. Внешним элементом является "весовая" версия списочного 
декодера Гурусвами-Судана для ОРС-кодов (см. алгоритм 17 из раздела 2), 
а внутренним — списочный переборный декодер кодов Адамара, кратко 
описанный в [3, с. 181]. Представим последний декодер в формализован- 
ном виде. 

Входными параметрами списочного переборного декодера кодов 


Адамара являются параметры (р”,т) -А-кода над полем /”, : мощность 
поля р, размерность кода 77 и упорядочение =21.--->=„ элементов про- 


Р 


странства А”,”. При декодировании на вход алгоритма подается слово 


У= ( У. ЕЙ ^ Декодер производит перебор всех кодовых слов и со- 
ставляет список У? их “весов по отношению к У”: 

АЛГОРИТМ 2: /* Вход: РТ; 2+1.---›2 „; У .Выход: \.*/ 
Шаг 0. Если массив С, кодовых слов (р” ,т) -А-кода пуст, то рассчитать 
его элементы: для всех 2,6 ЁР,, 1е {1;....р”} вычислить 
С. =\,„(2,); составить С. = (С.1>---- С, „), где \/„ - определено 
в 3.1; сохранить С для дальнейшего применения. 
Шаг 1. Для каждого 2,6 Р”, ЕЕ 4{1;.... р”} вычислить вес: 

и’ = тах{0;1- 4(>,С,,)Кр" - р”. 

Шаг 2. Составить вектор весов "’7 = (ил, 2-12 И’ т ) ивыдать У. 


Имея в наличии все необходимые алгоритмы, построим алгоритм 
списочного декодирования для КОРСА-кодов. Входными параметрами алго- 
ритма являются параметры КОРСА-кода: параметры полей Р и Т , длина 


Г и размерность А кода, упорядочения «*,>---> ® „„ и 2+1>--->2 „ эле- 


2 


ментов /^» и Ё’,” соответственно. При декодировании на вход алгорит- 


ма подается слово У = (у,.....У,)Е т Декодер производит поиск всех 
кодовых слов в пределах сферы, центром которой является у, радиусом — 
величина Е = (1-1/р”)(и- з/’р" (К /т- Т)). Выходом алгоритма яв- 


К 
ляется список всех информационных векторов Б(= Р’,), удовлетворяю- 


щих условию: 4(,„„.(Б), у) < Е ‚где У»... - кодирующее отображе- 
ние (”,А) -КОРСА-кода. Из [3, п.8.4.1] вытекает, что этот список содержит 
истинное сообщение. 
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АЛГОРИТМ 3: /* Вход: р.тьг.К :(3); 11 ›---> @ т Й И >--->- 7 ту 


У . Выход: список Б.*/ 
Шаг 0. Вычислить параметры: а = р”, & = К/т, к =г/а, 


в, = [/--0|, = а-ида- /&, - в’. 
Шаг 1. а) Разбить у = (>,,...у,) на блоки у’е Е?" : 
ода, р РЕЗ 


6) Для каждого /е {1;...;.} параметры Р, 1; 2.---->2„ и блок 


р 


у 


У подать на вход алгоритма 2; на выходе получить вектор весов 


и, = Сл >---> М7, ) . 


^ 


в) Составить вектор весов "’= ИР 5-5 МР ри + МР: МР, р 


сетку | | } | „| } и „| и | ‚| и 7 


Шаг 2. а) Параметры 94,7. Ао,Ё, , сетку У и вектор и’ подать на вход ал- 
горитма 1’. На выходе получить список {р,(х).....р,(х)}, где 
Р.Е Е» [1], РЕ, Де М 

6) Представить полиномы списка {р, (х),...,Р,(х)} в виде векторов: 
а. = Д„к(Рр, (ХЕ Е, ‚где Х,’‹ - определено в 3.1, #е {1;...;/}. Вы- 
дать список векторов {Ф.,...,Б,}, таких что 4(»,Б,) < Е, где 
Па 

4. Компьютерные модели списочных декодеров. 

4.1. Структурные схемы моделей списочных декодеров. На рис.1 и2 
изображены схемы декодеров и системы входов "весовой" версии списоч- 
ного декодера для ОРС-кодов и списочного декодера КОРСА-кодов. На вход 
декодера для ОРС-кодов (см.рис.1) подаются: мощность “Ч поля Галуа, 
длина Г и размерность А ОРС-кода, управляющий параметр 2, входное 
слово У и вектор весов у? помеченные на схеме как ВХ1-ВХб соответ- 


ственно. Декодирование происходит в блоках Б1-Б6. На вход декодера для 
КОРСА-кодов (рис.2) подаются параметры полей Галуа Р и 1, длина И 


и размерность А КОРСА-кода, и упорядочения & 1 и {2 элемен- 


тов ^„»” и Н,” соответственно и вектор У, помеченные на схеме как 
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ВХ1-ВХ7 соответственно. Декодирование происходит в блоках Б1-Б6. 
Рассмотрим подробнее работу блоков декодеров. 
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Декодер ОРС-кодов 





вх1 Вычисление параметров т, /, 

вх2 

вх3 

вх4 Б1 

вх5 — > Построение матрицы М 
вхб — № Б2 


Нахождение &л,л: 


Факторизация С в (Е[х]) [У] 
Су) =Х бу 


55 


1 Кох) 


Входы декодера: 
9 —_ ьч1 ГА —» вх4 


Рис.1. Структурная схема списочного декодера для ОРС-кодов. 


Схема входов 
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Декодер КОРСА-кодов 


вх1 > | Вычисление параметров К ”» 


вх2 
= 1,4, Е 
вх4 Б1 


Разбиение у на блоки у Ё”, 
Е Р 


вх7 —4 | уе Е 


Применение декодера кодов Адамара, 
получение: м -Е940,, )/а"-а" ие {1;... 
вр} | 


Построение о у=((@ 0), (0,0, )} 
вх5 »> и вектора весов \ = (№...) 
54 





список р(х) 


Формирование искомого списка 
екторов: ЕЕ”: ао>ь)<Е. 





+ список ВБ 


Входы декодера: 


ры > . 
„= 5 — № „к 
к > вх4 


Рис.2. Структурная схема списочного декодера 
для КОРСА-кодов. Схема входов 
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Работа блоков списочного декодера для ОРС-кодов состоит в следу- 
ЮЩЕем. Блок Б1 на вход получает 4,7’,А»# и вычисляет значения пара- 


Блок Б2 на вход получает /, 772, и строит матрицу однородной системы (2), 
обозначенную на схеме буквой М. Построение 1-й строки происходит путем 


вычисления коэффициентов (х,)”`” (у,)”” при неизвестной &..„ ( 


(хХьУрЕ {С У, ).....(х,.У,)Р) и присвоения их координатам строки. 
При этом порядок следования элементов строки несущественен, но в реа- 
лизации проще использовать естественный порядок, получаемый при рас- 
крытии двойной суммы. Блок БЗ на вход получает матрицу ЛМ и решает 
матричную систему ли = = о методом Гаусса, где = - вектор-столбец ко- 


эффициентов &,„,„ полинома (1), находит коэффициенты &„.,„. Блок 


54 получает на вход /,77,, коэффициенты @&,„_„ и строит полином (1). 


Блок Б5 реализует рекурсивную процедуру факторизации полинома (1) на 
основе алгоритма Рота-Руккенштейн. Процедура в процессе работы строит 
дерево коэффициентов полиномов, ветви которого образуют список эле- 
ментов, включающий искомый. Блок Бб кодирует полиномы списка выхода 
Б5, и формирует искомый список. 

Работа блоков списочного декодера для КОРСА-кодов состоит в следу- 
ющем. Блок Б1 на вход получает р,77,А, 7”, вычисляет значения парамет- 


ров декодера А, и, 1, = |= - У а, Е’. Блок Б2 на вход получает 7ъ и 


У , разбивает У на блоки >,, /Е {1:;...;,} для обработки алгоритмом 
2. Блок БЗ реализует алгоритм 2, вычисляя по полученным на вход блокам 


У, и упорядочению {2,}” веса \”’,,, где {Ее 41;...; р”\, 


1=1 
ДЕ \1;...;,}. Блок Б4, получая на вход веса %”,‚, и упорядочение 
{а,}”, ‚ формирует вектор весов у? и строит сетку У. Блок Б5, получая 
на вход У, \’, А, ,Ё применяет алгоритм 1’ и получает список полиномов 
Ко-1 
ВО а [х]. Блок Б6 получает на вход список полиномов и величину 


Е’, представляет полиномы в векторном виде и формирует искомый спи- 
сок. 
4.2. О программной реализации. Рассмотрим аспекты программной реа- 
лизации приведенных структурных схем. Вычисления в полях Галуа, век- 
торных пространствах и кольцах полиномов над полями Галуа реализованы 
на языке С-++ на базе динамической библиотеки \\ИпМТЕ-5_4_1 (см., напри- 
мер, [8]), включающей классы алгебраических структур и алгоритмов, 
необходимых для реализации моделей, таких как класс расширения поля 
Галуа, класс полиномов над полем и другие. Недостающие структуры и ал- 
горитмы получены в программной реализации в виде отдельных классов, 
например, кольцо полиномов с коэффициентами из кольца полиномов над 
полем Галуа. 

На основе полученной в предшествующей работе [5] реализации 
списочного декодера Судана для РС-кодов и указанных базовых компонен- 
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тов построены новые классы, реализующие рассмотренные структурные 
схемы декодеров. Для тестирования их работоспособности и постановки 
экспериментов получены реализации вспомогательных классов. 

Класс декодера Гурусвами-Судана для ОРС-кодов получен как на- 
следник класса списочного декодера Судана, включающий виртуальные 
процедуры создания улучшенной интерполяционной матрицы Гурусвами- 
Судана и обеспечения корректной работы с полями базового класса и клас- 
са-наследника с учетом специфики данного декодера. Класс декодера для 
КОРСА-кодов включает декодеры для ОРС-кодов и кодов Адамара как поля 
класса. 

Структурные схемы реализованы программно на основе библиотеки 
МЕС под следующие операционные системы: МИпаом!5 
95/98 /МТ/2000/ХР/\М а. 

Построенная программная реализация декодеров использована для 

проведения численных экспериментов в связи с применением списочного 
декодирования в схеме специального широковещательного шифрования 
[9], где имеет смысл использовать коды с относительно большим кодовым 
расстоянием. Так, например, при декодировании 140 слов (37,2)-ОРС-кода 
над полем А”, при числе ошибок в канале, не превышающем 70%, полу- 
чены списки объемом в одно кодовое слово, а в случае, когда число оши- 
бок составляло 70% — 81%, списки состояли из двух кодовых слов. Если 
число ошибок превышает 81%, то декодер не гарантирует правильное де- 
кодирование, так как при наших параметрах количество гарантируемо ис- 
правляемых ошибок равно [ х — \/"(А — 0—1] = 30 (см. раздел 2.). Деко- 
дирование 140 слов производилось программой в течение двенадцати се- 
кунд на компьютере с процессором мощностью 2,5 ГГц и ОЗУ объемом 512 
Мб. Из [9] вытекает, что рассмотренный пример ОРС-кода в схеме специ- 
ального широковещательного шифрования соответствует тиражу легально 
распространяемой продукции, равному 1369 экземпляров, а декодирование 
каждого слова гарантирует нахождение распространителей обнаруженного 
экземпляра контрафактной продукции. 
5. Заключение. Решены задачи разработки компьютерных моделей спи- 
сочных декодеров Гурусвами-Судана для ОРС-кодов и КОРСА-кодов с ис- 
пользованием эффективного алгоритма факторизации Рота-Руккенштейн. 
Для их применения в цифровых системах передачи данных к разработан- 
ной схеме можно добавить блок, реализующий выделение истинного сооб- 
щения из списка на выходе декодера [10]. На основе результатов, получен- 
ных в настоящей работе, возможно расширение компьютерной модели схе- 
мы специального широковещательного шифрования, построенной автором 
данной статьи в [9]. Отметим, что в настоящее время специалистами ве- 
дутся интенсивные теоретические исследования по оптимизации времени 
работы всех этапов списочного декодирования (см., например, [3], [11]). 
Разумеется, применение этих результатов в технической реализации деко- 
деров должно привести к улучшению их временных характеристик. 
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